


o 






X 



Phase transition and percolation 
in Gibbsian particle models 

Hans-Otto Georgii 

Mathematisches Institut der Universitat Miinchen 
Theresienstr. 39, D-80333 Miinchen, Germany. 



We discuss the interrelation between phase transitions in interacting lattice or contin- 
uum models, and the existence of infinite clusters in suitable random- graph models. In 
particular, we describe a random-geometric approach to the phase transition in the con- 
tinuum Ising model of two species of particles with soft or hard interspecies repulsion. 
We comment also on the related area-interaction process and on perfect simulation. 



, 1 Gibbs measures: general principles 

^^ , This section contains a brief introduction to the basic physical and stochastic ideas 

^^ I leading to the concept of Gibbs measures. The principal question is the following: 

Q>^ I Which kind of stochastic model is appropriate for the description of spatial 

^^ I random phenomena involving a very large number of components which are 

r^ ' coupled together by an interaction depending on their relative position? 

t^ ' 

To find an answer we will start with a spatially discrete situation; later we will proceed 

to the continuous case. Consider the phenomenon of ferromagnetism. A piece of ferro- 
magnetic material like iron or nickel can be imagined as consisting of many elementary 
magnets, the so-called spins, which are located at the sites of a crystal lattice and have 
a finite number of possible orientations (according to the symmetries of the crystal). 
The essential point is that these spins interact with each other in such a way that neigh- 
boring spins prefer to be aligned. This interaction is responsible for the phenomenon of 
spontaneous magnetization, meaning that at sufficiently low temperatures the system 
can choose between several distinct macrostates in which typically all spins have the 
same orientation. 

How can one find a mathematical model for such a ferromagnet? The first fact 
to observe is that the number of spins is very large. So, probabilistic experience with 
the law of large numbers suggests to approximate the large finite system by an infinite 
system in order to get clear-cut phenomena. This means that we should assume that the 
underlying crystal lattice is infinite. The simplest case to think of is the d-dimensional 
hypercubic lattice Z"'. (As the case d = 1 is rather trivial, we will always assume that 
d > 2.) On the other hand, to keep the model simple it is natural to assume that 
each spin has only finitely many possible orientations. In other words, the random 



spin ^j at lattice site i takes values in a finite state space S. The set of all possible 
spin configurations ^ = (Ci)jeZ'* is then the product space Q = S^ . This so-called 
configuration space is equipped with the Borel u-algebra JF for the natural product 
topology on 0,. Since the spins are random, we are interested in probability measures 
P on {Q,T). Such probability spaces are known as lattice systems. For any ^ G f2 and 
A C Z we write ^a = (^i)igA for the part of the configuration that occurs in A. By 
abuse of notation, we use the same symbol ^a for the projection from Q onto S . 

Which kind of probability measure on (J7, J^) can serve as a model of a ferromagnet? 
As we have seen above, the essential feature of a ferromagnet is the interaction between 
the spins. We are thus interested in probability measures P on fi for which the spin 
variables ^j, i G Z'^, are dependent. A natural way of describing dependencies is to 
prescribe certain conditional probabilities. This idea, which is familiar from Markov 
chains, turns out to be suitable also here. Since our parameter set Z admits no natural 
linear order, the conditional probabilities can, of course, not lead from a past to a future. 
Rather we prescribe the behavior of a finite set of spins when all other spins are fixed. 
In other words, we are interested in probability measures P on (0,^) having prescribed 
conditional probabilities 

GA(aiao (1) 

for a configuration ^a ^ S within a finite set A C Z"^ given a fixed configuration 
^A'^ £ S off A. In the following we write A CC Z when A is a finite subset of Z . The 
specific form of these conditional distributions does not matter at the moment. Two 
special cases are 

• the Markovian case: the conditional distribution (|l]) only depends on the value 
of the spins along the boundary dA = {i ^ A : |i — j| = 1 for some j € A} of A, 
i.e., 

GA{U\U^) = GA{a\^dA) (2) 

(with a slight abuse of notation); | • | stands for the Euclidean norm. 

• the Gibbsian case: the conditional distribution (Q) is defined in terms of a Hamil- 
ton function H^ by the Boltzmann-Gibbs formula 

GAieAiao = ^AiL ^M-HAm , (3) 

where Z^^^^ = Z]5'en:5'=^ off A ^'^v[~Ha[^')] is a normalizing constant. Physically, 
Ha{^) describes the energy excess of the total configuration ^ over the energy of 
the outer configuration ^a^- (Physicists will miss here the factor (3, the inverse 
temperature; we will assume that /3 is subsumed into Ha or, equivalently, that 
the units are chosen in such a way that /3 = 1.) 

In the following, Ga{ ■ |'^A'=) will be viewed as a probability measure on fi for which the 
configuration outside A is almost surely equal to ^a^- 

The above idea of prescribing conditional probabilities leads to the following concept 
introduced in the late 1960's independently by R.L. Dobrushin, and O.E. Lanford and 
D. Ruelle. 



Definition 1.1 A probability measure P on {0,,J^) is called a Gibbs measure, or DLR- 
state, for a fam,ily G = (GA)yvccZ'* of conditional probabilities (satisfying the natural 
consistency condition) if 

P{S,A occurs in A | ^a= occurs off A) = Ga(^a | ^a=) 

for P-almost all ^A'^ O'lT'd all A CC Z . 

If G is Gibbsian for a Hamiltonian i/ as in (|3|), eacli Gibbs measure can be interpreted as 
an equilibrium state for a physical system with state space 5 and interaction H. This 
is because the Boltzmann-Gibbs distribution maximizes the entropy when the mean 
energy is fixed; we will discuss this point in more detail in the continuum setting in 



Section [4.2| below. 

A general account of the theory of Gibbs measures can be found in the monograph 
Georgii (1988); here we will only present the principal ideas. In contrast to the situation 
for Markov chains, Gibbs measures do not exist automatically. However, in the present 
case of a finite state space S, Gibbs measures do exist whenever G is Markovian in the 
sense of (^), or almost Markovian in the sense that the conditional probabilities (||) are 
continuous functions of the outer configuration ^a^ • In this case one can show that any 
weak limit of Ga( • I^A^) for fixed ^ € il as A | Z"' is a Gibbs measure. 

The basic observation is that the Gibbs measures for a given consistent family G of 
conditional probabilities form a convex set Q. Therefore one is interested in its extremal 
points. These can be characterized as follows. 

Theorem 1.2 Let T = f]cr{^/^c : A CC Z ) the tail a-algebra, i.e., the a-algebra of 
all macroscopic events not depending on the values of any finite set of spins. Then the 
following statements hold: 

(a) A Gibbs measure P € Q is extremal in Q if and only if P is trivial on T , i.e., if 
and only if any tail measurable real function is P-almost surely constant. 

(b) Any two distinct extremal Gibbs measure are mutually singular on T . 

(c) Any non-extremal Gibbs measure is the barycenter of a unique probability weight 
on the set of extremal Gibbs measures. 

A proof can be found in Georgii (1988), Theorems (7.7) and (7.26). Statement (a) means 
that the extremal Gibbs measures are macroscopically deterministic: on the macroscopic 
level all randomness disappears, and an experimenter will get non-fluctuating measure- 
ments of macroscopic quantities like magnetization or energy per lattice site. Statement 
(b) asserts that distinct extremal Gibbs measures show different macroscopic behavior. 
So, they can be distinguished by looking at typical realizations of the spin configuration 
through macroscopic glasses. Finally, statement (c) implies that any realization which is 
typical for a non-extremal Gibbs measure is in fact typical for a suitable extremal Gibbs 
measure. In physical terms: any configuration which can be seen in nature is governed 
by an extremal Gibbs measure, and the non-extremal Gibbs measures can only be inter- 
preted in a Bayesian way as measures describing the uncertainty of the experimenter. 
These observations lead us to the following definition. 



Definition 1.3 Any extremal Gibbs measure is called a phase of the corresponding phys- 
ical system. If distinct phases exist, one says that a phase transition occurs. 

So, in terms of this definition the existence of phase transition is equivalent to the non- 
uniqueness of the Gibbs measure. In the hght of the preceding theorem, this corresponds 
to a "macroscopic ambivalence" of the system's behavior. We should add that not all 
critical phenomena in nature can be described in this way: even when the Gibbs measure 
is unique it may occur that it changes its qualitative behavior when some parameters 
of the interaction are changed. However, we will not discuss these possibilities here and 
stick to the definition above for definiteness. In this contribution we will ask: 

What are the driving forces giving rise to a phase transition? Is there any 
mechanism relating microscopic and macroscopic behavior of spins? 

As we will see, in a number of cases one can give the following geometric answer: 

One such mechanism is the formation of infinite clusters in suitably defined 
random graphs. Such infinite clusters serve as a link between the local and 
global behavior of spins, and make visible how the individual spins unite to 
form a specific collective behavior. 

In the next section we will discuss two lattice models for which this answer is correct. 
In Section 4 we will show that a similar answer can also be given in a continuum set-up. 
A useful technical tool is the stochastic comparison of probability measures. 

Suppose the state space 5 is a subset of R and thus linearly ordered. Then the 
configuration space fi has a natural partial order, and we can speak of increasing real 
functions. Let P, P' be two probability measures on Jl. We say that P is stochastically 
smaller than P' , and write P :< P', if / / dP < J f dP' for all local increasing functions 
(or, equivalently, for all measurable bounded increasing functions) / on fi. A sufficient 
condition for stochastic monotonicity is given in the proposition below. Although this 
condition refers to the case of finite products (for which stochastic monotonicity is 
similarly defined), it is also useful in the case of infinite product spaces. This is because 
(by the very definition) the relation ■< is preserved under weak limits. 

Proposition 1.4 (Holley's inequality) Let S be a finite subset ofH, A a finite index 
set, and P, P' two probability measures on the finite product space S^ giving positive 
weight to each element of S . Suppose the single-site conditional probabilities at any 
i £ A satisfy 

P{ ■ |Ca\{j} occurs off i) :< P'{ ■ \£,'\\u\ occurs off i) whenever ^ < ^' . 

Then P ^ P' . If this condition holds with P' = P then P has positive correlations in 
the sense that any two bounded increasing functions are positively correlated. 

For a proof (and a slight extension) we refer to Theorems 4.8 and 4.11 of Georgii, 
Haggstrom and Maes (1999). 



2 Phase transition and percolation: two lattice models 

To provide the necessary background for our results on continuum particle systems let 
us still stick to the lattice case. We will consider two classical models which allow an 
understanding of phase transition in random-geometric terms. Many further examples 
for the relation between random geometry and phase transition can be found in Georgii, 
Haggstrom and Maes (1999). 

Let us start recalling some basic facts on Bernoulli percolation on Z for d > 2. 
Consider Z*^ as a graph with vertex set Z'^ and edge set E{Z'^) = {e = {i,j} C Z'^ : 
\i — j\ = 1}. We fix two parameters < p.s,Pb ^ Ij the site and bond probabilities, and 
construct a random subgraph F = {X, E) of (Z'^, E{Z'^)) by setting 

X = {ieZ'^:ii = l}, E={ee E{X) : r/^ = 1} , 

where E{X) = {e E E{Z'^) : e C X} is the set of all edges between the sites of X, and ^i, 
i G Z"^, and r/e, e G E{Z'^), are independent Bernoulli variables satisfying P{ii = 1) = ps, 
P{r]e = 1) = Pb- This construction is called the Bernoulli mixed site-bond percolation 
model; setting pi, = 1 we obtain pure site percolation, and the case Ps = 1 corresponds 
to pure bond percolation. 

Let {0 <-> cxo} denote the event that T contains an infinite path starting from 0, and 

0(p„Pf,;Z'^) = Prob(Owoo) 

be its probability. By Kolmogorov's zero-one law, we have 6{ps,ph; Z'^) > if and only 
if r contains an infinite cluster with probability 1. In this case one says that percolation 
occurs. The following proposition asserts that this happens in a non-trivial region of 
the parameter square, which is separated by the so-called critical line from the region 
where all clusters of T are almost surely finite. The change of behavior at the critical 
line is the simplest example of a critical phenomenon. 

Proposition 2.1 The function 6{ps,ph;Z'^) is increasing in ps, pi, and d. Moreover, 
6{ps,Pb'-i Z'^) = when PsPb is small enough, while 0{ps,pi,; Z'^) > when d > 2 and PsPh 
is sufficiently close to 1. 



Sketch proof: The monotonicity in ps and pb follows from Proposition 1.4, and the one 
in d from an obvious embedding argument. To show that 6 = when PsPb is small, 
we note that the expected number of neighbors in P of a given lattice site is 2dpsPb- 
Comparison with a branching process thus shows that 9 = when 2dpsPi, < 1. 

Next, let d = 2 and suppose £ X but the cluster Cq of P containing is finite. 
Consider dextCo, the part of dCo belonging to the infinite component of Cq. For each 
site i € dextCo, either this site or all bonds leading from i to Cq do not belong to P. 
This occurs with probability at most 1 — PsPb- So, the probability that dextCo has a 
fixed location is at most (1 —psPbY with £ = #dextCQ. Counting all possibilities for this 
location one finds that 1 — 9 < 1 when 1 — PsPb is small enough. By the monotonicity 
in d, the same holds a fortiori in higher dimensions. □ 

The above proposition is all what we need here on Bernoulli percolation; an excellent 
source for a wealth of further rigorous results is the book of Grimmett (1999). 



We now ask for the role of percolation for Gibbs measures, and in particular for 
the existence of phase transitions. Of course, in contrast to the above this will involve 
dependent, i.e., non-Bernoulli percolation. We consider here two specific examples. In 
these examples, the family G of conditional probabilities will be Gibbsian for a nearest- 
neighbor interaction; this means that both (g) and (^) are valid. 

2.1 The Ising model 

This is by far the most famous model of Statistical Mechanics, named after E. Ising 
who studied this model in the early 1920s in his thesis suggested by W. Lenz. It is the 
simplest model of a ferromagnet in equilibrium. One assumes that the spins have only 
two possible orientations, and therefore defines S = {—1, 1}. The family G is defined 
by (ID with 

Ha{0 = J E %^?.}' (4) 

{J,j}nA7^0:|j-i| = l 

where J > is a coupling constant which is inversely proportional to the absolute tem- 
perature. This means that neighboring spins of different sign have to pay an energy cost 
J. There exist two configurations of minimal energy, so-called ground states, namely the 
configuration '+' which is identically equal to +1, and the configuration ' — ' identically 
equal to —1. The behavior of the model is governed by these two ground states. To 
see this we begin with some useful consequences of the ferromagnetic character of the 
interaction. First, it is intuitively obvious that the measures G^ = Ga(- |+) decrease 
stochastically when A increases (since then the effect of the + boundary decreases). 



This follows easily from Holley's inequality, Proposition 1.4. Since the local increasing 
functions are a convergence determining class, it follows that the weak infinite-volume 
limit P"*" = lim^|2'* G*^ exists. Likewise, the weak limit P~ = lim^|2<^ G'a(- |— ) exists 
(and by symmetry is the image of P"*" under simultaneous spin fiip). These limits are 
Gibbs measures and invariant under translations. Holley's inequality also implies that 
they are stochastically maximal resp. minimal in G, and in particular extremal. This 
gives us the following criterion for phase transition in the Ising model. 

Proposition 2.2 For the Ising model on Z with any coupling constant J > we have 
#G > I if and only if P^ 7^ P+ if and only if J ^0 dP^ > 0. 

The last equivalence follows from the relation P^ :< P^ , the translation invariance of 
these Gibbs measures, and the spin-fiip symmetry. A detailed proof of the proposition 
and the previous statements can be found in Section 4.3 of Georgii, Haggstrom and 
Maes (1999). 

How can we use this criterion? This is where random geometry enters the scenery. 
The key is the following geometric construction tracing back to Fortuin and Kasteleyn 
(1972) and in this form to Edwards and Sokal (1988). It is called the random- cluster 
representation of the Ising model. 

Let 8^ = {E C E{Z'^) : E D E{k'')} be the set of all edge configurations in Z"^ 
which include all edges outside A, and define a probability measure 0a on 8'^ by setting 

0a(P) = ^a' 2'=(^)p#^\^(^^) {l-p)#Ei^')\E when E D P(A-), (5) 



where p = 1 — e~'^ , k{E) is the number of clusters of the graph (Z'^,£'), and Z^ is 
a normaUzing constant. (/>a is called the random- cluster distribution in A with wired 
boundary condition. This measure turns out to be related to G\ = G\(- |+). It will 
be convenient to identify a configuration ^ G J7 with the pair (X+,X~), where X+ and 
X~ are the sets of all lattice sites i for which ^j = +1 resp. —1. 

Proposition 2.3 For any hypercube A in Z there exists the following correspondence 
between the the Gibbs distribution G^ for the Ising model and the random- cluster dis- 
tribution (/)A in ([3j. 

( G~j^ -^ (pA ) Take a spin configuration ^ = {X^,X~) G Q, with distribution G^, and 
define an edge configuration E G S'^ as follows: Independently for all e € i?(Z°') let 
e G -E with probability 



1 



PA(e) 



- e-^ if e C A:+ or e C X' 


, and e n A 7^ 


1 if e C A^ 




otherwise, 





and e ^ E otherwise. Then E has distribution <j)\. 

(0A '^ G^ ) Pick an edge configuration E G £^ according to (/>a, and define a spin 
configuration ^ = {X~^,X~) G 0, as follows: For each finite cluster C of {Z'^,E) let 
G C X^ or G C X^ according to independent flips of a fair coin; the unique infinite 
cluster of {Zi'^,E) containing A^ is included into X^ . Then ^ has distribution G^. 

Proof. A joint description of a spin configuration ^ G $1 with distribution G^ and 
an edge configuration E G £j^ with distribution (p\ can be obtained as follows. Any 
edge e G E{Zi ) is independently included into E with probability p = 1 — e~ resp. 1 
according to whether e H A 7^ or not; each spin in A is equal to +1 or —1 according 
to independent flips of a fair coin; the spins off A are set equal to +1. The measure 
P thus described is then conditioned on the event A that no spins of different sign are 
connected by an edge. Relative to P{ ■ \A), ^ has distribution G^ and E has distribution 
(/)A. This is because exp[— //a('?)] is equal to the conditional P-probability of A given 
^, and 2^^^^'^ is proportional to the conditional P-probability of A given E. Now, it 
is easy to see that the two constructions in the proposition simply correspond to the 
conditional distributions of E given ^ resp. of S, given E relative to P{ ■ \A). □ 

Intuitively, the edges in the random-cluster representation indicate which pairs of spins 
"realize" their interaction, in that they decide to take the same orientation to avoid 
the dealignment costs. On the one hand, this representation is the basis of an efficient 
simulation procedure, the algorithm of Swendsen-Wang (1987), which together with its 
continuous counterpart will be discussed at the end of Section [4.4| . On the other hand, 
it is the key for a geometric approach to the phenomenon of phase transition, as we will 
now show. 

The construction (p\ -^ Gj implies that, for G A, the conditional expectation of 
.^0 given £^ is 1 when is connected to dA by edges in E, and otherwise. Hence 

J^odPl = MO^dA). (6) 



By Proposition 1.4, the measures (j)\ decrease stochastically when A increases, so that 



the infinite-volume random-cluster distribution (p = lim^^2;<^ '/'A exists. Letting A | Z"^ 



we 



in (|6|) we thus find that / ^o dP^ = 0(0 ^^ oo) . Combining this with Proposition 2/2 
obtain the first statement of the following theorem, the equivalence of percolation and 
phase transition. This gives us detailed information on the existence of phase transition. 

Theorem 2.4 Consider the Ising model on Z with Hamiltonian ^ for any coupling 
constant J > 0. Then ^Q > 1 if and only if (p{0 <-> oo) > 0. Consequently, there exists 
a coupling threshold < Jc < oo (corresponding to a critical inverse temperature) such 
that ^Q = 1 when J < Jc and ^Q > 1 when J > Jc- 

Sketch proof. It only remains to show the second statement. This follows from Holley's 



inequality, Proposition L4. First, this inequality implies that 0a is stochastically in- 
creasing in the parameter p = 1 — e~^ . Hence 0(0 ^^ oo) is an increasing function of p. 
Moreover, one finds that is stochastically dominated by the Bernoulli bond percola- 
tion measure, whence 0(0 <->• oo) = when p is so small that 0{l,p; Z ) = 0. Finally, 
is stochastically larger than the Bernoulli bond percolation measure with parameter 
p = p/{p + 2(1 — p)). Hence 0(0 <-> oo) > when p is so large that 6{l,p; Z'^) > 0. 
Details of this computation can be found in Section 6 of Georgii, Haggstrom and Maes 
(1999), which deals in fact with the extension of these results to the Potts model in 
which each spin has q different values. □ 

One may ask whether the connection between percolation and phase transition can be 
seen more directly in the behavior of spins. The following corollary gives an answer to 
this question. Let {0 < — > oo} denote the event that belongs to an infinite cluster of 
the graph {X~^ , E{X^)) induced by the set of plus-spins. 

Corollary 2.5 For the Ising model on Z with arbitrary coupling constant J > we 
have P^{0 < — > oo) > whenever ^Q > 1. The converse holds only when d = 2. 



The first part follows readily from the construction in Proposition 2.3 which shows 
that P+(0 < — > oo) > 0(0 (^ oo). For its second part see Georgii, Haggstrom and 
Maes (1999). Pursuing the idea of plus-percolation further one can obtain the following 
result independently obtained in the late 1970s by Aizenman and Higuchi on the basis 
of previous work of L. Russo; a simpler proof has recently been given by Georgii and 
Higuchi (1999). 

Theorem 2.6 For the Ising model on T? with J > Jc, there exist no other phases than 
P+ andP-. 

A celebrated result of Dobrushin asserts that in three or more dimensions there exist 
non-translation invariant Gibbs measures which look like P"*" in one half-space and like 
P^ in the other half-space. 

2.2 The Widom-Rowlinson lattice gas 

The Widom-Rowlinson lattice gas is a discrete analog of a continuum model to be 
considered in Section 4. It describes the random configurations of particles of two 



different types, plus or minus, which can only sit at the sites of the lattice 7/" . Multiple 
occupations are excluded. So, at each site i of the lattice there are three possibilities: 
either i is occupied by a plus-particle, or by a minus-particle, or i is empty. The 
configuration space is thus $7 = S"^ with S = {—1,0,1}. The basic assumption is 
that there is a hard-core repulsion between plus- and minus-particles, which means that 
particles of distinct type are not allowed to sit next to each other. In addition, there 
exists a chemical "activity" z > which governs the overall-density of particles. The 
Hamiltonian thus takes the form 

Ha{0 = E U{C^,CJ) - log 2 E l^^l ' (7) 

{i,i}nA^0:|i-j|=i «eA 

where C^(CiiG') ~ "-^ '^^ dCj = ~1) ^^^ U{^i,^j) = otherwise. The associated family G 
of conditional probabilities is again given by (|^). Just as in the Ising model, for z > 1 
there exist two distinguished configurations of minimal energy, namely the constant 
configurations '-I-' and '— ' for which all sites are occupied by particles of the same type. 
Moreover, one can again apply Holley's inequality to show that the Gibbs distributions 
G\ = Ga{ -1+) and G^ = Ga( -I—) converge to translation invariant limits P^,P~ G G 
which are stochastically maximal resp. minimal in G, and therefore extremal. This 
implies that Proposition |2.2| holds verbatim also in the present case. 

Is there also a geometric representation of the model, just as for the Ising model? 
The answer is yes, with interesting analogies and differences. There exists again a 
random-cluster distribution with an appearance very similar to (^), but this involves 
site percolation rather than bond percolation. Namely, consider the probability measure 
V'A on the set X^ = {Y C Z'^ : Y D A^)} which is given by 

MY) = ^A ^ 2'^(^) p#^\^^ (1 - p)#2'\^ for FDA- (8) 

here p = j§r7 ^(^) is the number of clusters of the graph (Y,E{Y)), and Z\ is a nor- 
malizing constant. V'A is called the site random- cluster distribution in A with parameter 
p and wired boundary condition. We will again identify a configuration ^ € il with a 
pair (X^ ,X^), where X^ and X~ are the sets of all lattice points i such that ^j = +1 
resp. —1; thus X^ U X~ is the set of occupied sites, and ^ = on its complement. 
Here is the random- cluster representation of the Widom-Rowlinson lattice gas which is 
analogous to Proposition |2.3|. 



Proposition 2.7 For any hypercube A in TJ^ there exists the following correspondence 
between the the Gibbs distribution G^ for the Widom-Rowlinson model and the site 
random- cluster distribution ip\ in (^. 

(GJ -^ ^A ) For a random spin configuration ^ = {X^,X^) E 0, with distribution 
G^, the random set Y = X^ U X^ has distribution ip\. 

{ip\ -^ G^ ) Pick a random set Y E X^ according to f/'A; (md define a spin con- 
figuration ^ = {X^,X^) E il with X^ U X^ = Y as follows: For each finite cluster 
G of {Y,E{Y)) let G C X+ or C C X^ according to independent flips of a fair coin; 
the unique infinite cluster of (Y,E(Y)) containing A^ is included into X^ . Then ^ has 
distribution Gt ■ 



The random-cluster representation of the Widom-RowUnson model is simpler than that 
of the Ising model because the randomness involves only the sites of the lattice. (This is 
a consequence of the hard-core interaction; in the case of a soft repulsion the situation 
would be different, as we will see in the continuum setting in Section 4.) On the other 
hand, there is a serious drawback of the site random-cluster distribution ijj\: it does 



not satisfy the conditions of Proposition 1.4 for positive correlations. This is because 



the conditional probabilities in (|T^) below are not increasing in Y . So, we still have a 
counterpart to @, viz. 

l^odGt = ^PA{0^dA), (9) 

but we do not know if the measures V'A are stochastically increasing in z and stochas- 
tically decreasing in A. So we obtain a somewhat weaker theorem. 

Theorem 2.8 Consider the Widom-Rowlinson model on Z , d > 2, defined by l^ with 
activity z > 0. Then #Q > I if and only i/lim^i^'' V'aIO ^^ dA) > 0. In particular, we 
have ^Q = 1 when z is sufficiently small, and ^Q > 1 when z is large enough. 

Sketch proof. The first statement follows immediately from (M) and the analog of Propo- 



sition 2.2. To prove the second statement we note that Va has single-site conditional 



probabilities of the form 

^PAiY 3i\Y\ {i}) = p/\p +{l-p) 2''(^'^)-i] , (10) 

where «;(i, Y) is the number of clusters of y \ {i} that intersect a neighbor of i. Since 
< k(z, Y) < 2d, it follows from Proposition |LJ that ■0A is stochastically dominated by 
the site-Bernoulli measure with parameter p* = p/{p + (1 — p)2~^), and dominates the 
site-Bernoulli measure with parameter p* = p/{p+ (1 — p)2^'^~^). Combining this with 
Proposition 2.1 we thus find that ^Q = 1 when z is so small that 9{p* , 1; Z"^) = 0, and 
#^ > 1 when z is so large that 9{p^, 1; Z'') > 0. □ 



Let us note that, in contrast to Theorem p.4| , the preceding result does not extend to the 
case when there are more than two different types of particles; this is related to the lack 
of stochastic monotonicity in this model. However, due to its simpler random-cluster 
representation the Widom-Rowlinson model has one advantage over the Ising model, in 



that it satisfies a much stronger counterpart to Corollary 2.5. In analogy to the notation 
there we write {0 < — > oo} for the event that the origin belongs to an infinite cluster of 
plus-particles. 

Corollary 2.9 For the Widom-Rowlinson lattice gas on Z for arbitrary dimension 
d >2 and with any activity z > 0, we have ^G > ^ if and only if P^(0 < *• oo) > 0. 



Sketch proof. The construction in Proposition 2.7 readily implies that 

^a(0 ^ dA) = GXiO ^^ dA) . 

Combining this with (^) and letting A | Z one obtains the result. □ 

The above equivalence of phase transition and percolation even holds when Z is re- 
placed by an arbitrary graph. 
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As noticed before Theorem |2.8| , we have no stochastic monotonicity in the activity z, 
and therefore no activity threshold for the existence of a phase transition. We are thus 
led to ask if, at least, the particle density is an increasing function of z. This can be 
deduced from general thermodynamic principles relying on convexity of thermodynamic 
functions rather than stochastic monotonicity. This will be described in Section 4.2 in 
the continuum set-up. 

3 Continuum percolation 

In the rest of this contribution we will show that quite a lot of the preceding results and 
techniques carry over to models of point particles in Euclidean space. In this section 
we deal with a simple model of continuum percolation. Roughly speaking, this model 
consists of Poisson points which are connected by Bernoulli edges. To be precise, let X 
denote the set of all locally finite subsets X of R . X is the set of all point configurations 
in R , and is equipped with the usual u-algebra generated by the counting variables 
X -^ i^X\ for bounded Borel sets A C R , here we use the abbreviation Aa = A n A. 
Next, let £ be the set of all locally finite subsets of £'(R'^) = {{x, y} C R*^ : x / y}. £ is 
the set of all possible edge configurations and is equipped with an analogous a-algebra. 
For A € -^ let E{X) = {e G E(R'^) : e C A} the set of aU possible edges between the 
points of A, and £x = {E & £ : E C E{X)} the set of edge configurations between the 
points of A. We construct a random graph T = (A, E) in R'^ as follows. 

• Pick a random point configuration A € Af according to the Poisson point process 
vr^ on R with intensity z > 0. 

• For given X £ X, pick a random edge configuration E € £x according to the 
Bernoulli measure /J^ on £x for which the events {E B e}, e = {x,y} S E{X), 
are independent with probability /^^(-E B e) = p{x — y); here p : R"^ —^ [0, 1] is a 
given even measurable function. 

The distribution of our random graph F is thus determined by the probability measure 

P''P{dX,dE)=7r'{dX)fiP^{dE) (11) 

on <Y X £". It is called the Poisson random- edge model, or Poisson random- connection 
model, and has been introduced and studied first by M. Penrose (1991); a detailed 
account of its properties is given in Meester and Roy (1996). 

A special case of particular interest is when p{x — y) is equal to 1 when \x — y\ < 2r 
for some r > 0, and otherwise. This means that any two points x and y are connected 
by an edge if and only if the balls Br{x) and Br{x) with radius r and center x resp. y 
overlap. The connectivity properties of the corresponding Poisson random-edge model 
are thus the same as those of the random set H = [j^^x -^rix), for random A with 
distribution vr^. This special case is therefore called the Boolean model, or the Poisson 
blob model. 

Returning to the general case, we consider the percolation probability of a typical 
point. Writing x ^-> oo when x belongs to an infinite cluster of F = (A, E), this is given 
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by the expression 

»(,, ,., R^) = j *'""^^^-"°°' P'^'idX. iE) (12) 

for an arbitrary bounded box A with volume |A|. By translation invariance, 9{z,p;'Rr) 
does not depend on A. In fact, in terms of the Palm measure tt^ of tt^ and the associated 
measure P^'P{dX,dE) = ■JT^{dX) fixidE) we can write 

e(z,p;R^) = P"'P(0^oo) . 

The following result of M. Penrose (1991) is the continuum analog of Proposition p.l| . 

Theorem 3.1 9{z,p; R ) is an increasing function of the intensity z and the edge prob- 
ability function p{-) . Moreover, 6{z,p;'R,'^) = when z J p{x) dx is sufficiently small, 
while 9{z,p;H^) > when z J p{x) dx is large enough. 

Sketch proof: The monotonicity follows from an obvious stochastic comparison argu- 
ment. Since z J p{x) dx is the expected number of edges emanating from a given point, 
a branching argument shows that 6[z,p;'R,'^) = when z J p{x) dx < 1. It remains to 
show that 9{z,p;Il ) > when z J p{x) dx is large enough. By scaling we can assume 
that J p{x) dx = 1. For simplicity we will in fact suppose that p is bounded away from 
in a neighborhood of the origin, i.e., p{x — y) > S > whenever |a; — y| < 2r. (The 
following is a special case of an argument of Georgii and Haggstrom (1996).) 

We divide the space R'^ into cubic cells A(i), i G Z*^, with diameter at most r. We 
also pick a sufficiently large number n and introduce the following two concepts. 



• 



Call a cell A(i) good if it contains at least n points which form a connected set 
relative to the edges of F in between them. This event does not depend on the 
configurations in all other cells and has probability at least 



TT^iNi > n) [1 - (n - 1)(1 - 52)n-2] = p^ 



here, Ni is the random number of points in cell A(i), and the second term in the 
square bracket is an estimate for the probability that one of the n points is not 
connected to the first point by a sequence of two edges. The essential fact is that 
Ps is arbitrarily close to 1 when n and z are large enough. 

• Call two adjacent cells A{i),A{j) linked if there exists an edge from some point 
in A(i) to some point in A(j). Conditionally on the event that A(i) and A(j) are 
good, this has probability at least 1 — (1 — 8)"' = pi,, which is also close to 1 when 
n is large enough. 

Now the point is the following: whenever there exists an infinite cluster of linked good 
cells (i.e., an infinite cluster in the countable graph with vertices at the good cells and 
with edges between pairs of linked cells) then there exists an infinite cluster in the 
original Poisson random-edge model. Hence 9{z,p; R'^) > - 9{ps,Pb', Z'^), where v is the 
cell volume and 9{ps,Pb','Zi'^) is as in Proposition |2.1| . Hence 9{z,p;'R'^) > when z is 
large enough. D 
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How can one extend a percolation result as above from the Poisson case to point pro- 
cesses with spatial dependencies? Just as in the lattice gas, one can take advantage of 
stochastic comparison techniques. To this end we need a continuum analog of Holley's 
theorem. 

A simple point process P on a bounded Borel subset A of R'^ (i-e., a probability 
measure on A'a = {X € A:" : X C A}) is said to have Papangelou (conditional) intensity 
7 : A X X\ — > [0, oo[ if P satisfies the identity 

/ P{dX) Y. fix, X \ {x}) = fdx f P{dX) 7(xlA) fix, X) (13) 

for any measurable function / : A x X\ -^ [0, oo[. (This is a non-stationary ana- 
log of the Georgii-Nguyen-Zessin equality discussed in the contribution of D. Stoyan 
to this volume.) This equation roughly means that 7(x|X) dx is proportional to the 
conditional probability for the existence of a particle in an infinitesimal volume dx 
when the remaining configuration is X. Formally, it is not difficult to see that (p!^) 
is equivalent to the statement that P is absolutely continuous with respect to the 
intensity-1 Poisson point process tta = vrj^ in A with Radon-Nikodym density g sat- 
isfying giX U {x}) = 7(x|X)g((X), see e.g. Georgii and Kiineth (1997). In particular, 
the Poisson process vr^ of intensity z > on A has Papangelou intensity jix\X) = z. 

Proposition 3.2 (Holley-Preston inequality) Let A C R'^ be a bounded Borel set 
and P, P' two probability measures on X\ with Papangelou intensities 7 resp. 7'. Suppose 
jix\X) < j'ix\X') whenever X C X' and x ^ X' \X . Then P ■< P' . If this condition 
holds with P' = P then P has positive correlations. 

Of course, the stochastic partial order P ^ P' is defined by means of the inclusion 
relation on X\. Under additional technical assumptions the preceding proposition was 
first derived by Preston in 1975; in the present form it is due to Georgii and Kiineth 
(1997). In the next section we will see how this result can be used to establish perco- 
lation in certain continuum random-cluster models, and thereby the existence of phase 
transitions in certain continuum particle systems. 

4 The continuum Ising model 

The continuum Ising model is a model of point particles in R of two different types, 
plus and minus. Rather than of particles of different types, one may also think of 
particles with a ferromagnetic spin with two possible orientations. The latter would be 
suitable for modelling ferrofluids such as the Au-Co alloy, which have recently found 
some physical attention. Much of what follows can also be extended to systems with 
more than two types, but we stick here to the simplest case. A configuration of particles 
is then described by a pair ^ = (A+,X^), where X^ and X^ are the configurations of 
plus- resp. minus-particles. The configuration space is thus Q, = X'^ . 

We assume that the particles interact via a repulsive interspecies pair potential 
of finite range, which is given by an even measurable function J : R -^ [0, 00] of 
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bounded support. The Hamiltonian in a bounded Borel set A C R*^ of a configuration 
^ = {X~^,X~) is thus given by 

Ha{0= E Jix-y). (14) 

xeX+ ,y£X- : {x,y}nA^li 

In view of its analogy to @ this model is called the continuum Ising model. Setting 
J{x — y) = oo when |x — y| < 2r and J{x — y) = otherwise, we obtain the classical 
Widom-Rowlinson model (1970) with a hard-core interspecies repulsion (which in spatial 
statistics is occasionally referred to as the penetrable spheres mixture model). Of course, 
this case corresponds to the Widom-Rowlinson lattice gas considered above. Here 
we make only the much weaker assumption that J is bounded away from zero on a 
neighborhood of the origin. That is, there exist constants 6,r > such that 

J(x — y) > 5 when |a; — y| < 2r. (15) 

An interesting generalization of the Hamiltonian (O) can be obtained by adding 
an interaction term which is independent of the types of the particles. In a ferrofluid 
model this would mean that in addition to the ferromagnetic interaction of particle 
spins there is also a molecular interaction which is spin-independent. Such an extension 
is considered in Georgii and Haggstrom (1996). 

Given the Hamiltonian (|l^), the associated Gibbs distribution in A with activity 
z > and boundary condition (^j^c = (X^c,Xj^c) G -^Ic is defined by the formula 

GAidUlU^) = ^A|J^, expi-H^iO] ^AidXt) TTlidXX) (16) 

which is completely analogous to (^. The corresponding set G = Q{z) of Gibbs measures 
is then defined as in Definition p^ (with A running through the bounded Borel sets in 
R"' instead of the finite subsets of Z'^). 

In general, the existence of Gibbs measures in continuum models is not easy to 
establish. In the present case, however, it is simple: Thinking of X];^ as the space of 
configurations on two disjoint copies A"*" and A~ of A, we see that Ga{ ■ I^A^) has the 
Papangelou intensity 

T(:.|X+,X-) = p»P[-p»«--'<^-';] ""f |<.^ (17) 

[ z exp[-J2yex+J{^-y)\ iixeA J 

Proposition ^^ therefore implies that G'a(-|^ac) ^ vr^ x vr^. Standard compactness 
theorems for point processes now show that for each ^ G Af^ the sequence Ga{- I^a^) 
has an accumulation point P as A | R , and it is easy to see that P £ Q. 

4.1 Uniqueness and phase transition 

We will now show that the Gibbs measure is unique when z is small, whereas a phase 
transition occurs when z is large enough. Both results rely on percolation techniques. As 
in the Widom-Rowlinson lattice gas, it remains open whether there is a sharp activity 
threshold separating intervals of uniqueness and non- uniqueness. 
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Proposition 4.1 For the continuum Ising model we have 4j^Q{z) = 1 when z is suffi- 
ciently small. 

Sketch proof: Let P,P' € G{z). We show that P = P' when z is small enough. Let R 
be the range of J, i.e., J{x) = when |a;| < R, and divide R'^ into cubic cells A(i), 
i S Z*^, of linear size R. Let p* be the Bernoulli site percolation threshold of the graph 
with vertex set Z*^ and edges between all points having distance 1 in the max-norm. 
Consider the Poisson measure Q^ = ir^ x vr^ on the configuration space 17 = X^. 

Let ^,^' be two independent realizations of Q^, and suppose z is so small that 
Q^ X Q^{Ni + Nl > 1) < p*, where Ni and A^^' are the numbers of particles (plus or 
minus) in ^ resp. ^'. Then for any finite union A of cells we have Q^ x Q^{A <^— > oo) = 0, 
where {A <^— > oo} denotes the event that a cell in A belongs to an infinite connected set 



of cells A(i) containing at least one particle in either ^ or ^'. Proposition 3.2 together 



with ( |T7|) imply that P x P' < Q^ x Q^ . Hence P x P'{A ^^-^ oo) = 0. In other words, 
given two independent realizations ^ and ^' of P and P' there exists a random corridor 
of width R around A which is completely free of particles. In particular, this means 
that ^ and ^' coincide on this random corridor. By a spatial strong Markov property of 
Gibbs measures, it follows that P and P' coincide on the u-algebra of events in A. As 
A can be chosen arbitrarily large, this proves the proposition. □ 

After this result on the absence of phase transition (following from the absence of some 
kind of percolation) we turn to the existence of phase transition. This will follow 
from the existence of percolation in a suitable random-cluster model. In analogy to 



Propositions ^^ and ^J, we will derive a random-cluster representation of the Gibbs 
distribution 

G+ = y'^Xc(dy^+)GA(-|y+,0) (18) 

with a Poisson boundary condition of plus-particles and no minus-particle off A. Its 
random-cluster counterpart is the following probability measure XA on JY x £ describing 
random graphs (1", E) in R : 

XAidY, dE) = Z-|V^^ 2^(^'^) T:\dY) f^p^{dE) . (19) 

In the above, k{Y,E) is the number of clusters of the graph {Y,E), Z/^^Yj^c = 
J 2^^^'^' 7T^{dY/^) normalizes the conditional probability of XA given Y/^c = y n A'^ 
(so that Y\c still has the Poisson distribution vr^c), and fiy i^ ^^^ probability measure 
on £ for which the edges e = {x, y} C Y are drawn independently with probability 
p{x — y) = 1 — e''^^^~'^' if e ^ Ya.c, and probability 1 otherwise. The probability mea- 
sure XA in (0) is called the continuum random- cluster distribution in A with connection 
probability function p and wired boundary condition. Note that (in contrast to (S) and 
(§)) this distribution describes random configurations of both points and edges. In 
the Widom-Rowlinson case of a hard-core interspecies repulsion the randomness of the 
edges disappears, and XA describes a dependent Boolean percolation model which is the 
direct continuum analog of (^. The random- cluster representation of the continuum 
Ising model now reads as follows. 
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^-J(x-y) if e c x+ or eCX- 


, and e n A 7^ 


1 if e C A^ 




otherwise. 





Proposition 4.2 For any bounded box A in R"^ there is the following correspondence be- 
tween the Gibbs distribution Gj^ in pq j for the continuum, Ising model and the random- 
cluster distribution XA in ^1^). 

( G J -^^ Xa) Take a particle configuration ^ = {X^,X~) G Q with distribution 
Gjy and define a random graph {Y,E) £ X x £ as follows: Let Y = X^ U X^ , and 
independently for all e = {x, y} € E{Y) let e £ E with probability 



PA{e) 



Then {Y,E) has distribution XA- 

( XA '^ Gj^ ) Pick a random graph (Y, E) £ X x £ according to XA ■ Define a particle 
configuration S, = {X^,X~) G Q with X^ U X^ = Y as follows: For each finite cluster 
C of (y, E) let C C X^ or C d X^ according to independent flips of a fair coin; 
the unique infinite cluster of {Y-, E) containing Y\c is included into X^ . Then £, has 
distribution G\. 

Just as in the lattice case, the random-cluster representation above gives the following 
key identity: for any finite box A C A, 

= J#{xeY^:x^ Y^c} XA{dY, dE) ; (20) 

in the above, the notation x <-> Ya^ means that x is connected to a point of Y\c in the 
graph (y, £'). In other words, the difference between the mean number of plus- and 
minus-particles in A corresponds to the percolation probability in xa- How can one 
check that the latter is positive for large z? The idea is again a stochastic comparison. 
Let v\ = xa(- X £) the point marginal of XA- Then XA{dY, dE) = v\{dY) (j)\y{dE) 
with an obvious analog </>A,y of (|5[) . An application of Proposition |1.4| shows that <J)a,y 
is stochastically larger than the Bernoulli edge measure ^y ^^ which edges are drawn 
independently between points x,y £ Y with probability p{x — y) = (1 — e~ )/[(l — 
e~ ) + 2e~ ] when jx — y| < 2r, and with probability otherwise. Here, 6 and r are as 
in assumption ([T5|). Moreover, z^a has the Papangelou intensity 

7(x|y) = z /2'=(^uW'-)d<^A,yuw/ f ^'^"^''^ dc^A 



Yu{x}/ / ^ ^ ' ' acpAy . 

To get a lower estimate for 'y{x\Y) one has to compare the effect on the number of 
clusters in (Y, E) when a particle at x and corresponding edges are added. In principle, 
this procedure could connect a large number of distinct clusters lying close to x, so that 
k{Y U {x}, • ) was much smaller than /c(y, • ). However, one can show that this occurs 
only with small probability, so that ^{x\Y) > az for some a > 0. By Proposition 



3.2| , we can conclude that XA is stochastically larger than the Poisson random-edge 



measure P^^^'P defined in (0). The right-hand side of ( po|) is therefore not smaller than 
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9{az,p;'R'^). Finally, since Gj^ ^ vr^ x vr^ by (0), the Gibbs distributions Gj have a 
cluster point P~^ G Q{z) satisfying 

|[#X+ -#X^] P+{dX+,dX-) > e{az,p;-R'') . 

By spatial averaging one can achieve that P^ is in addition translation invariant. To- 



gether with Theorem 3.1 this leads to the following theorem. 

Theorem 4.3 For the continuum Ising model on YC^, d > 2, with Hamiltonian (il 



and sufficiently large activity z there exist two translation invariant Gibbs measures P^ 
and P~ having a majority of plus- resp. minus-particles and related to each other by the 
plus-minus interchange. 

This result is due to Georgii and Haggstrom (1996). In the special case of the Widom- 
Rowlinson model it has been derived independently in the same way by Chayes, Chayes, 
and Kotecky (1995). The first proof of phase transition in the Widom-Rowlinson model 
was found by Ruelle in 1971, and for a soft but strong repulsion by Lebowitz and Lieb 
in 1972. Gruber and Griffiths (1986) used a direct comparison with the lattice Ising 
model in the case of a species-independent background hard core. 

As a matter of fact, one can make further use of stochastic monotonicity. (In contrast 
to the preceding theorem, this only works in the present case of two particle types.) 
Introduce a partial order '<' on $7 = X'^ by writing 

{X+,X-) < {Y+,Y-) when X+ C Y+ and X' D ¥'. (21) 



A straightforward extension of Proposition 3.2 then shows that the measures Gj^ in 
( p!8|) decrease stochastically relative to this order when A increases. (This can be also 
deduced from the couplings obtained by perfect simulation, see Section ^^ below.) It 



follows that P"*" is in fact the limit of these measures, and is in particular translation 
invariant. Moreover, one can see that P~^ is stochastically maximal in Q in this order. 
This gives us the following counterpart to Corollary |2.9|. 



Corollary 4.4 For the continuum Ising model with any activity z > 0, a phase transi- 
tion occurs if and only if 



f P+{dX+,dX-)fiP^+{0 ^^ oo) > 



here P'^ is the Palm measure of P^ , and the relation < — > oo means that the origin 
belongs to an infinite cluster in the graph with vertex set X~^ and random edges drawn 
according to the probability function p = I — e~^ . 

It is not known whether P+ and P^ are the only extremal elements of Q{z) when 
d = 2, as it is the case in the lattice Ising model. However, using a technique known in 
physics as the Mermin-Wagner theorem one can show the following. 

Theorem 4.5 If J is twice continuously differentiable then each P € Q{z) is translation 
invariant. 

A proof can be found in Georgii (1999). The existence of non-translation invariant 
Gibbs measures in dimensions d > 3 is an open problem. 
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4.2 Thermodynamic aspects 

Although we were able to take some advantage of stochastic comparison techniques in 



the continuum Ising model, the use of Proposition |3.2| is much more limited than that 
of its lattice analog. The reason is that its condition requires some kind of attractivity, 
which is in conflict with stability (preventing the existence of infinitely many particles 
in a bounded region). This implies that a continuum Gibbs distribution G\ with a pair 



interaction cannot satisfy the conditions of Theorem 3.2 with P = P' = G\, which 
would imply that G\ has positive correlations and is stochastically increasing with the 
activity z. Fortunately, this gap can be closed to some extent by the use of classical 
convexity techniques of Statistical Mechanics. These will allow us to conclude that at 
least the particle density of Gibbs measures is an increasing function of the activity z. 

It should be noted that these ideas are standard in Statistical Physics; they are 
included here because they might be less known among spatial statisticians, and because 
we need to check that the general principles really work in the model at hand. One 
should also note that this technique does not depend on the specific features of the 
model; in particular, it applies also to the gas of hard balls discussed in H. Lowen's 
contribution to this volume. 

Let us begin recalling the thermodynamic justification of Gibbs measures. Let A C 
R be a finite box. For any translation invariant probability measure P on $7 = Af^ 
consider the entropy per volume 

s(P)= hm \k\-^S{PK). 

[A|^oo 

Here we write Pa for the restriction of P to Xf^, the set of particle configurations in A, 

and 

^ f - / log / dPfy if Pa < vrj^ X ■k\ with density /, 

1 — oo otherwise 

for the entropy of Pa relative two the two-species Poisson process vr^ x tt\ on Af| with 
intensity 1. The notation |A| -^ oo means that A runs through a specified increasing 
sequence of cubic boxes with integer sidelength. The existence of s(P) is a multidimen- 
sional version of Shannon's theorem; see Georgii (1988) for the lattice case to which the 
present case can be reduced by identifying O with {Xq)'^ for a unit cube C. 
Next consider the interaction energy per volume 

u(P)= /'p(dX+,dX-)[l{oex+} E Ji^) + h^ex-} E ^^^ (22) 

xex- X&X+ 

defined in terms of the Palm measure P of P. u(P) can also be defined as a per- volume 
limit, cf. Georgii (1994), Section 3. Also, consider the particle density 

q{p) = p{n) = \A\-' J[#x+ + #xx] dPK 

of P; by translation invariance the last term does not depend on A. The term 
—q{P) \ogz is then equal to the chemical energy per volume. 
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Finally, consider the pressure 

p(z) = - min [u(P) - g{P) log z - s(P)l ; (23) 

the minimum extends over all translation invariant probability measures P on 0,. The 
large deviation techniques of Georgii (1994) show that p{z) = lim|^|_^oQ l^l"*^ log-^AI^Ac 
for each ^ G $7. (The paper Georgii (1994) deals only with particles of a single type and 
superstable interaction, but the extension to the present case is straightforward because 
J is nonnegative and has finite range.) The variational principle for Gibbs measures 
then reads as follows. 

Theorem 4.6 Let P be a translation invariant probability measure on Q = X'^ . Then 
P G Q{z) if and only if u{P) — q{P) log z — s{P), the free energy per volume, is equal to 
its minimum —p{z). 

The "only if" part can be derived along the lines of Georgii (1994) and Proposition 7.7 
of Georgii (1995). The "if" part follows from the analogous lattice result (see Section 
15.4 of Georgii (1988)) by the identification of i7 and (X^) mentioned above. 

What does the theorem tell us about the particle densities of Gibbs measures? Let 
us look at the pressure p{z). First, it follows straight from the definition ( p3| ) that p(z) 
is a convex function of logz. In other words, the function p(t) = p(e*) is convex. Next, 



inserting P = vr^ x Jg into the right-hand side of (23) we see that p > — oo, and that 
the slope of p at t tends to infinity as t ^ co. 

Now, suppose P G G{z) is translation invariant. The variational principle above 
then implies that the function t —>■ {t — log z)p{P) + p(z) is a tangent to p at log z. For, 
on the one hand we have 

u(P) - q{P) logz- s(P) = -p(z) < oo 

and thus u(P) — s(P) < cxd, and on the other hand 

u(P) - q{P) t - s(P) > -p(t) for all t. 

Inserting the former identity into the last inequality we get the result. As a consequence, 
the particle density g{P) lies in the interval between the left and right derivative of p 
at log 2. By convexity, these derivatives are increasing and almost everywhere identical. 
In fact, they are strictly increasing. For, if they were constant on some non-empty open 
interval / then for each tQ £ I and P € Q{e*°) the function t -^ Q{P)t — p(t) would 
be constant on /, and thus by the variational principle P G ^(e*) for all t £ I. This 
is impossible because the conditional Gibbs distributions depend non-trivially on the 
activity. We thus arrive at the following conclusion. 

Corollary 4.7 Let < z < z' and P € G^z), P' £ Q{z') be translation invariant. Then 
q{P) < q{P'), and g{P) ^ oo as z —s- oo. 

In the present two-species model it is natural to consider also the case when each 
particle species has its own activity, i.e., the plus-particles have activity z~^ and the 
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minus-particles have activity z~ . It then follows in the same way that the pressure 
p(z^,z^) is a strictly convex function of the pair (logz^,logz~), and therefore that 
the density of plus-particles is a strictly increasing function of 2^, and the density of 
minus-particles is a strictly increasing function of z~ ; these densities tend to infinity as 



z^ resp. z~ tends to infinity. Moreover, Theorem 4.3 implies that p(z"'',z~) has a kink 
at [z^z) when z is large enough; this means that the convex function t -^ p(ze*,ze~^) 
is not differ entiable at t = 0. 

4.3 Projection on plus-particles 

As the continuum Ising model is a two-species model, it is natural to ask what kind of 
system appears if we forget all minus-particles and only retain the plus-particles. The 
answer is that their distribution is again Gibbsian for a suitable Hamiltonian. This 
holds also in the case of different activities z~^ and z~ of plus- and minus-particles, 
which is the natural context here. To check this, take any box A C R and let A D A 
be so large that the distance of A from A'^ exceeds the range R of J. Integrating over 
X^ in (^) and conditioning on X'^l^.^ one finds that the conditional distribution of X'^ 
for given ^^\^ under Ga{ ■ |Ca=) does not depend on A and ^ac and has the Gibbsian 
form 

G'A(dA+|X+.) = Z-i^^^ exp[-z- H^{X+)] tt^ (dXX) (24) 

with the Hamiltonian 

Ha{X+) = / (1 - exp [ - ^ J(x - y)]) dy . (25) 

•'^ XGX+ 

Thus, writing Q(z~^,z') for the set of all continuum-Ising Gibbs measures on il = X^ 
with Hamiltonian (^) and activities z~^ and z^ , and Q{z^,z^) for the set of all Gibbs 
measures on JY with conditional distributions (p3), one obtains the following corollary. 



Corollary 4.8 Let P E G{z^ , z ) and P be the distribution of the configuration of plus - 
particles. Then P € G^z^, z^). In particular, ^Q{z,z) > 1 when z is large enough. 



The last statement follows from Theorem 4.3. 



In the Widom-Rowlinson case when J = 00 In . |<2j.}, the relationship between 
Q{z^,z~) and Q{z^,z^) has already been observed in the original paper by Widom 
and Rowlinson (1970). The plus-Hamiltonian ( |25| ) then takes the simple form 



^A(A+) = lAn U B2r{x)\ 
X&X+ 

that is, Ha[X^) is the volume in A of the Boolean model with radius 2r induced 
by X^ . In this case, the model with distribution G /^{dX'^\X'^c) was reinvented by 
Baddeley and van Lieshout (1995). Having the two-dimensional case in mind, they 
coined the suggestive term area-interaction process. From here one can go one step 
further to Hamiltonians which use not only the volume but also the other Minkowski 
functionals. This has been initiated by Likos et al. (1995) and Mecke (1996); see also 
Mecke's contribution to this volume. 
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One particularly nice feature of the area-interaction process and its generalization 
5D is that it seems to be the only known (non-Poisson) model to which Proposition 



3.2 can be applied for establishing positive correlations of increasing functions. This 
attractiveness property makes the model quite attractive for statistical modelling. (By 
way of contrast, repulsive point systems can be modelled quite easily, for example by a 
nonnegative pair interaction.) 

However, some caution is necessary due to the phase transition when z^ = z~ = z 
is large: The typical configurations of Ga( -10) for a large finite window A then can be 
typical for either phase, P^ or P^, and thus can have different particle densities. Due 
to finite size effects, this phenomenon already appears when z+ and z" are sufficiently 
close to each other. So, the spatial statistician should be aware of such an instability of 
observations and should examine whether this is realistic or not in the situation to be 
modelled. 



Finally, one can use Proposition 3^ to show that the Gibbs measures P G G{z'^ ^ z^) 
are stochastically increasing in z+ and decreasing in z^ . In particular, the density of 
plus-particles for any P E Q{z'^ , z^) increases when z+ increases or z~ decreases, as can 



also be seen using the partial order (^11). The monotonicity results in the last paragraph 



of Section 4.2 thus follow also from stochastic comparison techniques, but Corollary 4.7 



cannot be derived in this way. 

4.4 Simulation 

There are various reasons for performing Monte-Carlo simulations of physical or sta- 
tistical systems, as discussed in a number of other contributions to this volume. In the 
present context, the primary reason is to sharpen the intuition on the system's behav- 
ior, so that one can see which properties can be expected to hold. This can lead to 
conjectures which then hopefully can be checked rigorously. 

Here we will show briefly how one can obtain simulation pictures of the contin- 
uum Ising model. We start with a continuum Gibbs sampler which is suggested by 



Proposition 4.2; in the Widom-Rowlinson case it has been proposed by Haggstrom, van 
Lieshout and M0ller (1997). 

Consider a fixed window A C R'' and the Gibbs distribution G^^^ = Ga(- |0, 0) 
for the continuum Ising model in A with activity z and free (i.e., empty) boundary 
condition off A. We define a random map F : X\ — > <Ya by the following algorithm: 

• take an input configuration X G X\, 

• select a Poisson configuration Y € X\ with distribution vr^, 

• define a random edge configuration E C {{x, y} : x G X, y € 1"} by independently 
drawing an edge from x ^ X to y ^Y with probability p{x — y) = 1 — e~'^^^~y\ 



set F{X) = ly £Y ■.{x,y} ^E for all x G x\. 



That is, F{X) is a random thinning of vr^ obtained by removing all points which 
are connected to X by a random edge. Its distribution is nothing other than the 
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Poisson point process tt^W on A with inhoniogeneous intensity measure pi^{dy) = 
z 1a (y) exp [— X]a;GX "^ iv ~ ^)]} '^V ■ (^^ course, this could also be achieved by setting 
F{X) = {y G Y : Uy < exp [— ^^gj,^ J(y — x)]} for independent [7(0, l)-variables Uy, 
y ^Y. Although this was simpler in the case of the MCMC below, it would considerably 



increase the running time of the perfect algorithm of Theorem 4.10, as will be explained 



there.) Now the point is that tTa'i -5^+ is the conditional distribution of X~ given X~^ 
relative to the Gibbs distribution G^'^'^, and similarly with + and — interchanged. So, 
if ^ = (X+,X~) has distribution G^*^^ and F^,F~ are independent realizations of F 
then {F~^{X~),F~' o F^{X^)) has again distribution G^"'^. This observation gives rise 
to the following Markov chain Monte Carlo algorithm (MCMC). 

Proposition 4.9 Let F^,F~ , n >0, be independent realizations of F, and Xq S X\ 
any initial configuration. Define recursively 

X, = F,-iX+), X+ = F+iX-_,), X- = F-{X+) for n > 1. 

Then the distribution of {X^ ,X~) converges to G^^^^ in total variation norm at a geo- 
metric rate. 

Proof. It suffices to observe that F = with probability 5 = e~^ ' ' . This shows that 
for any two configurations AT, X' € Af^. and any A C <^a 

|Prob(F(A:) G A) - Prob(F(X') G A)| < 1 - 5 . 

So, if one looks at the process (A'+,X~) for two different starting configurations then 
each application of F reduces the total variation distance by a factor of 1 — (5. □ 

A nice property of the random mapping F is its monotonicity: if X C X' then 
F{X) D F{X') almost surely. This allows to modify the preceding algorithm to obtain 
perfect simulation in the spirit of Propp and Wilson, as described in the contribution 
of E. Thonnes to this volume. According to (p^, G^^^'^ is stochastically dominated by 
independent Poisson processes of plus- and minus-particles. So one can use the idea of 
dominated perfect simulation in her terminology. We describe here only the algorithm 
and refer to her contribution for more details. 

Roughly speaking, the perfect algorithm consists of repeated simultaneous runs of 
the preceding MCMC, starting from two particular initial conditions at some time N}^ < 
until time 0. The two initial conditions are chosen extremal relative to the ordering 
(^T]), namely with no initial plus-particle (the minimal case), and with a Poisson crowd 
of intensity z of plus-particles (which is maximal by stochastic domination). Since the 
same realizations of F are used in both cases, the two parallel MCMC's have a positive 
chance of coalescing during the time interval from Nk to 0. If this occurs, one stops. 
Otherwise one performs a further run which starts at some time A'^^.+i < Nk- 

Theorem 4.10 Let F^,F^ , n < 0, be independent realizations of F, and (A'^fc)fc>i a 
strictly decreasing sequence of negative run starting times. For each run indexed by 
k>l let 

^k = F+oFZ,oF^,o...oF^^^,oF+^^,oF^^ 
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be the random mapping corresponding to the MCMC of Proposition ^ for the time 
interval from Nk to 0, and consider the processes 

Xt^n = ^'^ (0) ' ^fctmax = $fc O F+ (0) . 

Then there exists a smallest (random) K < oo such that Xj^ ^^^ = Xj^ ^^^, and the 
random particle configuration ^k = (-^Kmin'-^o ("'^/rmin)) ^^^ distribution G^*^*^. 

Since the random mapping F can be simulated by simple standard procedures, the 
implementation of the preceding algorithm is quite easy; a Macintosh application can be 
found at tittp : //www . mathematik . uni-muenchen . de/~georgii/CIsing . html| . The 



main task is to store the random edge configuration E in each application of F during 
a time interval {N^,- ■ ■ ^N^-i — 1} for use in the later runs (which should be done 
in a file on the hard disk when z\K\ is large). As a matter of fact, once the set E is 
determined one can forget the positions of the particles of Y and only keep their indices. 
In this sense, E contains all essential information of the mapping F. As a consequence, 
knowing E one needs almost no time to apply the same realization of F in later runs. 
This is not the case for the alternative definition of F mentioned above. 

However, there are some difficulties coming from the phase transition of the model. 
Running the perfect algorithm for small z is fine and raises no problem. But if z is large 
then the algorithm requires a considerably longer time to terminate. This is because 
for each run k the distribution of ^fc,max = (-'^A^max'-^o~(^Atmax)) "^i^^ ^e close to P+ 
and thus show a large crowd of plus-particles giving the minus-particles only a minimal 
chance to spread out. Likewise, the distribution of ^fc.min = (^A^'min' -^o~(^fc'min)) ^^^^ 
be close to P~, so that the minus-particles are in the great majority. The bottleneck 
between these two types of configurations is so small that K will typically be much 
too large for practical purposes, at least for windows A of satisfactory size. In order 
to reduce this difficulty, one should not simulate the Gibbs distribution G^^'^ with free 
boundary condition (as we have done above for simplicity), as this distributes most of 
its mass on two quite opposed events. Rather one should simulate one of the phases, 
say P^ . For a finite window, this can be achieved by imposing a random boundary 



condition of Poisson plus-particles outside A as in (18). Such a boundary condition 
helps -'^^jnax ^^^ "'^fc'min Q^ite a lot to coalesce within reasonable time. (If one is willing 
to accept long running times, one should impose periodic boundary conditions to reduce 
the finite-size effects.) 

This electronic version does not contain the figures, 
which can be found on the website mentioned above. 

The pictures shown are obtained in this way. The underlying interaction potential 
is J(x) = 3(1 — |x|)^ for \x\ < 1, J(x) = otherwise. The size of the window is 
20 X 20. Outside of the window there is an invisible boundary condition of white Poisson 
particles. The activities are z = 4.0 (top) and z = 4.5 (bottom), and the corresponding 
coalescence times are —Nk = 300 resp. 400. 

In the subcritical case z = 4, the particles in the bulk do not feel the boundary 
condition: there is no dominance of white over black. In the supercritical case z = 4.5 
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however, the influence of the white boundary condition is strong enough to dominate 
the whole window, and the phase transition becomes manifest. This is nicely illustrated 



by the random-cluster representation, which according to Proposition 4.2 is obtained 
from the point configuration {X~^,X~) by adding random edges within X~^ and X^ 
separately. Here one sees that in the subcritical case the influence of the plus-boundary 
condition is only felt by a the particles near the boundary, while in the supercritical 
case the global behavior is dominated by a macroscopic cluster reaching from the bound- 
ary far into the interior of A. This visualizes the equivalence of phase transition and 



percolation derived in Corollary 4.4. 



To conclude we mention two other algorithms. First, there is another perfect algo- 
rithm using a rejection scheme due to Fill, which has been studied in detail by Thonnes 
(1999) in the Widom-Rowlinson case; its extension to the present case is straightfor- 
ward. A further possibility, which is particularly useful in the supercritical case, is to 
use a continuum analog of the Swendsen-Wang algorithm, Swendsen and Wang (1987). 
In its classical version for the lattice Ising model, this algorithm consists in alternating 



applications of the two procedures in Proposition 2.3. In its continuum version, one can 
again alternate between the two procedures of Proposition 4.2, but one has to combine 
this with applications of the random mapping F in order to obtain a resampling of par- 
ticle positions. Unfortunately, this algorithm does not seem to admit a perfect version 
because of its lack of monotonicity, but it has the advantage of working also for the 
many-species extension of the model. 
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